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, Abstract 

The rare decay B ^ Xsjj is studied in the Standard Model (SM) and in two different 
^ ■ versions (Model I and Model II) of the Two Higgs Doublet Model (2HDM). In the SM the 

, branching ratio into hard photons is about 1 x 10~^ and can be appreciably different in 

the 2IIDM. We also introduce a forward-backward asymmetry which gives an additional 
handle to discriminate different models. 
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5 ■ 1 Introduction 

p • Intense experimental effort is being directed to B-physics. Several new facilities are on the 
I horizon. There are, of course, the e~^e~ based symmetric and asymmetric B-factories at Cornell, 
KEK and SLAC, in addition to LEP II. Progress is also being made in the hadronic environment 
^ I at HERA-B and there are some plans for TEV-B and far into the future, LHC-B. Indeed it is 
useful to recall that in the case of the kaon some of the branching ratios currently being measured 
have reached the 10~^° level. It should then be clear that in B-physics too, the experimental 
activity is likely to continue to flourish well beyond the presently attainable branching ratios 
of about 10"^. 

Amongst the rare decays the flavor changing decays of the B meson are of special interest, in 
particular those driven by the electroweak (EW) penguins, due to their relative cleanliness and 
to their sensitivity to new physics [Q . In this category, the two decays that have dominated the 
discussion for over a decade are b —>■ sj and b —>■ sl^l^ . It is useful to recall that the inclusive 
branching ratio for b ^ s'y has been measured at CLEO to be 

Br{B X,7)^^P = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10"^ (1) 
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and a Next-to-Leading Order calculation of the same now exists The result 0]: Br{B 

^slY^ = (3.28 ± 0.33) X 10~^ indicates agreement of the SM prediction with the CLEO mea- 
surement within 2a. As far as 6 ^ is concerned, the non-resonant part of -B ^ Xsfi~^fi~ 
has been predicted to give 0: Br{B Xgji^ ji^)]^^ = (5.7 ± 0.9) x 10"^, while the existing 
upper bound on this decay mode from D0 is now being updated to 

Br{B < 3.2 x 10"^ , (2) 

less than one order of magnitude away from the theoretical prediction. 

Our focus in this paper will be on a related mode: b 577 which is of the same order in the 
EW couplings as 6 sl^l~ and is also of considerable interest. We will extend earlier works 
on this process and study it in the SM and in one of its most popular extensions, namely the 
2HDM. The calculation of b —>■ 577 in the SM at the lowest order in the EW interactions and 
without QCD corrections gives a branching ratio of about 10^'', therefore in the ballpark of the 
rare B decays that should become accessible at the future B-meson facilities. Furthermore, we 
will introduce a forward-backward asymmetry which should be less sensitive to QCD corrections 
than the branching ratio and should also constitute a useful probe of the theory. 

2 The b S77 decay in the Standard Model 

The decay B ^^770 can be studied to a very good approximation in terms of the quark level 
decay b 577 Q. The total amplitude for the quark level process can be written as [|ToHl^ 

A{b .77) = -t^^e,{h)e,{k2)u{ps)n{pb) , (3) 
V 27r 

where ae is the electromagnetic fine structure constant, Gp the Fermi coupling constant, pb 
and Ps denote the momenta of the incoming and outgoing quarks and e^(A;i) and e,y(/c2) are the 
polarization vectors of the two photons. The tensor T^'^ is derived from the calculation of the 
Feynman diagrams of Fig. |l], when we sum over the three possible flavors of quarks that run in 
the penguin loop 

T'^-= J2 XiTr = Xu{Tr -TD + MTr -TD ■ (4) 

i=u,c,t 

^The inclusive calculation presented here is applicable to B^, Bd and Bs with the corresponding multiparticle 
state {Xs) having the flavors of su, sd and ss respectively. Moreover our calculations can be readily adapted 
to the case of 5 ^ ^77 with obvious changes. 
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Figure 1: Examples of IPI and IPR contributions to the b —>■ s'-ff process. The dashed hne 
represents a W (SM) or a charged scalar (2HDM). 



Here Aj = Vif,V*. is the CKM factor multiplying the loop diagrams with an internal i-quark and 
we have used the unitarity of the CKM matrix to obtain the last result. Even though <^ At, 
for our purpose, as we will explain later, the first term cannot be neglected. 

As we see from Fig. 1, both one-particle reducible (IPR) and one-particle irreducible (IPI) 
diagrams contribute to the process. We will follow the notation of Ref. |jl2| and write the total 
amplitude for the i flavor as 



j,lPI + j,lPR 



(5) 



The IPI contribution T.^^^pj is then written as 



7;Tpi = 9^3(^0/'^'^ , 



(6) 



where Zi = 2ki ■ k2/mf and 5^{zi) is a function defined by ||T3| 



53(2;^) = 1 + — f — log [1 - ZiUil - u) 
Zi Jo u 



(7) 



The tensor /^"^ is given by0 



.kM:ko 



te^'^'-jaliiki - k2)^ + zf^(e'^5.«^. _ e^^'^^k^)^^jL 

Ki ■ K2 



(8) 



where we use the notation 'Jl,r = (1 =F75)/2. On the other hand, the IPR contribution has the 
form 



^iPR = -^lF2{xi)W''''{m,^R + m,7i) 



(9) 



^We want to emphasize here that both Refs. [|iO[ and deal with the decay Bs 77. In that case one 
only needs to keep the first term in the /^"^ tensor. This is not true for the inclusive B Xs"fj decay. 
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where F2{xi) for Xi = mf/M^ is the form factor of the magnetic dipole operator of the 657 
vertex [|l^], given by 



F2(Xi) 



[-AQ + 205xi-312x^i + 175x^i-22xf + {36x^i-5Axf)\og{xi))] , (10) 



36(1 

while the W^" tensor is defined by the expression 



Pb 



+ 2 



Ps -h Pb- k2 
1 1 



a{fi, ki 



P's 



Pb 



Ps ■ k2 Pb - ki 



a{u, k2)a{^, ki) + 



Ps -ki Pb- ki 

1 1 



cr(z/, k2) 



Ps - ki Pb - k2 



The rate r{b S77) can therefore be decomposed as the sum of a pure IPR, a pure IPI 
and an interference IPR-lPI term, i.e. 

T{b -> S77) = FipR + Pipi + PipR.ipi . (12) 

In order to obtain the total rate into hard photons we have to place suitable kinematical cuts and 
perform one dimensional integration numerically. We have also checked our results integrating 
over the phase space with a Montecarlo event generator. 

Indeed the total rate and the relevance of each different contribution (1PR,1PI) depends 
appreciably on the kinematical cuts imposed. In order to isolate two hard photons, we demand 
that their energy is not too small and that they are not too collinear to each other and to the 
outgoing s quark. We thus require the energy of each photon to be larger than 100 MeV and 
that Si = {ps + kiY, $2 = {ki + ^2)^ and S3 = {ps + ^2)^ satisfy 



Si> cm^ , S2 > cm^, 



S3 > crrif, 



(13) 



We take c = 0.01 and c = 0.02 to study the dependence on the cuts. Note that the resulting two 
photon invariant mass squared is at least one order of magnitude bigger than m^. Furthermore 
all the angles are taken to be > 20°. 

Our results are summarized in Table |I|. These results are for pure EW penguins and do not 
include QCD corrections [1^. For rrit ~ 175 GeV, as in the b —>■ s'j case, QCD corrections are 
not likely to change the predictions of Table |l| dramatically. The branching ratios are calculated 
by normalizing to the semileptonic branching ratio, for which we have used the experimental 
value Br{B X^ez/e) ~ 0.11. We did not include QCD corrections in the semileptonic rate 
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Br(6^ 


S77) X 10 










Total 


1-Particle 


1-Particle 


Interference 


^FB 






Reducible 


Irreducible 








c=0.01 


SM 




1.60 


1.30 


0.23 


0.08 


0.69 




tan P 












2HDM 


0.5 


0.23-14.67 


0.003-14.70 


0.23 


-(0.26-0.003) 


0.37-0.77 


Model I 


1 


0.23-1.26 


0.01-0.96 


0.23 


-0.007-0.07 


0.39-0.67 




10 


1.57-1.59 


1.26-1.29 


0.23 


0.07 


0.68 


2HDM 


0.5 


2.19-16.67 


1.87-16.17 


0.23 


0.09-0.27 


0.70-0.75 


Model II 


1 


2.07-9.65 


1.75-9.21 


0.23 


0.09-0.21 


0.70-0.74 




1 Cl 


2.03-7.76 


1.71-7.34 


0.23 


0.09-0.18 


0.70-0.74 






c=0.02 


SM 




1.33 


1.02 


0.23 


0.08 


0.64 




tan (3 












2HDM 


0.5 


0.23-11.57 


0.003-11.60 


0.23 


-(0.26-0.003) 


0.37-0.73 


Model I 


1 


0.23-1.05 


0.009-0.76 


0.23 


-0.007-0.07 


0.38-0.62 




10 


1.30-1.33 


0.99-1.01 


0.23 


0.07 


0.64 


2HDM 


0.5 


1.80-13.26 


1.48-12.76 


0.23 


0.09-0.27 


0.66-0.71 


Model II 


1 


1.70-7.71 


1.38-7.23 


0.23 


0.09-0.21 


0.65-0.71 




10 


1.67-6.21 


1.35-5.79 


0.23 


0.089-0.18 


0.65-0.70 



Table 1: Values of branching ratio and forward-backward asymmetry obtained in the SM and 
in the 2HDM's (Model I and II), for two different values of the cut. In the 2HDM case we give 
a range of values corresponding to the variation of between 100 GeV and 1 TeV. 

either; for our purpose QCD effects are completely marginal here. For the numerical results 
we have used: rrit = 175 GeV, rric = 1.5 GeV, m„ = 4.5 McV and = 4.5 GcV. However, we 
have checked the dependence of the results on rric and m^, allowing them to vary in the ranges 
1.2 < rric < 1.8 GeV and 4.5 < < 100 MeV. We find that the dependence on m„ is totally 
negligible while the branching ratio varies up to 15 — 20% within this very conservative range 
of nic- For we have made the following distinction: in the calculation of the quark decay 
amplitude we have taken the current quark mass, = 150 MeV, while in the integration over 
the phase space, in order to respect the kinematics of the decay as closely as we can, we have 
used the K meson mass and taken approximately = 450 MeV. Again this assumption does 
not affect the results by more than 8% in any case. As far as the CKM matrix elements are 
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concerned, we have used: Vcs — Vtb — 1, — 0.22, \Vts\ — \Vcb\ — 0.04 and |V„fe/Kf<| — 0.08. 

In Table we also give results separately due to IPR, IPI and due to their interference. As 
we can see, the role played by the IPR and IPI contributions varies a little with the different 
cuts imposed. With mild cuts, the IPR contribution is dominant, due to the divergent behavior 
of the amplitude for S2 ^ 0. This divergent behavior is caused by the presence of both infrared 
{E-y — > 0) and coUinear (cos 6*^^ — > 0) divergencies, which, for sure, would cancel in a complete 
analytical calculation when virtual photon corrections at the same order in are also included. 
Experimentally, if one of the two photons in 6 — S77 decay is emitted with very low energy 
and/or at a very small angle with respect to the other photon, it becomes difficult to distinguish 
the event from a 6 — S7 decay. Therefore we exclude this region of the phase space by imposing 
the cuts mentioned earlier. This tends to suppress the IPR contribution compared to the IPI 
one. 

0.045 
0.04 
0.035 

<c^ 0.03 
T3 

^ 0.025 
T3 

^ 0.02 
0.015 
0.01 
0.005 
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A 

S2 

Figure 2: The normalized distribution l/F • dV /ds2 versus S2 = S2/ml in the SM for two values 
of the cut parameter: c = 0.01 (solid) and c = 0.02 (dashed). Also, (52) = 0.25,0.28 for 
c = 0.01orc = 0.02 respectively. 




Moreover, the analytical expressions of the IPR and IPI terms in Eqs. (^-(pA|) are such 
that the light quarks mainly contribute to the IPI part and the heavy quarks to the IPR one. 



This is due to the analytical behavior of the coefficient functions 6s{zi) and F2{xi) (see Ref. [12 



for a plot of these functions). The physical reason for this is that, when the top quark runs 
in the loop of the penguin diagram for b s'-ff (both in the IPR and in the IPI diagrams), 
the amplitude for the decay is very well approximated by the amplitude for the emission of a 
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brehmstrahlung photon off the s quark leg. Indeed by this reasoning, the authors of Ref. 



were able to get the right expression for the top contribution using Low's theorem |]16|. Thus 
the higher the kinematical cuts imposed to reduce the s quark brehmstrahlung events, the 
higher the relevance of the light quark contribution (i.e. IPI) to the total amplitude. For this 
reason we retained in Eq. (^) both the Xt and the A„ terms, in order to have the light quark 
dependence fully under control. 
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Figure 3: The normalized distribution 1 /T ■ dT / d{cos versus cos^^^ in the SM for two values 
of the cut parameter: c = 0.01 (solid) and c = 0.02 (dashed). Also, {cos9^^) = —0.70,-0.68 
for c = 0.01 or c = 0.02 respectively. 

In Fig. 1^ and |^ we present the results for two different distributions: dV /ds2 {s2 = S2/ml) 
and dT / d{cos 9^^) , where 9^^ is the angle between the two photons, for the two sets of cut^. As 
we can see the kinematics of the process is such that the s quark tends to be emitted with rather 
high energy, compensated by the harder of the two photons, while the less energetic photon 
tends to go in the direction of the s quark. This topology is typical of a brehmstrahlung event 
of the s quark and in fact it gets enhanced when we choose milder cuts on the energies and on 
the angles. 

Although the b 577 process is higher order in ae and rarer than b S7, it does allow 

one to introduce an additional feature that may be a useful probe of the dynamics of the decay, 

namely a forward-backward asymmetry, defined as 

•^We observe that the behavior of 1/T ■ dT/ds2 around §2 — 0.45 (i.e. S2 — 9GeV^) corresponds, for ttIc — 1.5 
GeV, to S2 — 4m^, i.e. to the physical threshold for an on-shell cc pair. 
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A 



FB 



r(cos0,^ > 0) - r(cos0,^ < 0) 



(14) 



r(cos Os^t > 0) + r(cos < 0) ' 

where Os-^ is the angle between the s quark and the softer photon. Indeed, due to the identical 
nature of the two photons, this is the only non-trivial angle we can think of to study a forward- 
backward asymmetry. We recall that a forward-backward asymmetry has also been found to 
be useful in the study of the h — > sl^l~ decay Jl^. Our results are summarized in Table |l], 
for two different values of the cuts. In passing, we must remark that we expect A-pB to be less 
affected by QCD corrections than the branching ratio. 



3 The h 577 decay in the Two Higgs Doublet Model 

We want now to consider the h — >■ 577 decay in the context of a 2HDM with no flavor changing 
neutral currents allowed at the tree level, i.e. Model I and Model II, in which the up and down 
type quarks couple respectively to the same or to two different Higgs doublets. 




100 200 300 400 500 600 700 800 900 1000 

Mh (GeV) 



Figure 4: The ratio Rbt as a function of the charged scalar mass Mh, for different values of 
tan/3 in Model I (solid: tan/5 = 0.5; long-dashed: tan/3 = 1; short-dashed: tan/3 = 10) and 
Model II (dotted: tan/3 = 0.5; long-dash-dot: tan/3 = 1; short-dash-dot: tan/3 = 10). The 
curves correspond to the cut parameter c = 0.01. 

The physical charged scalar field contributes to the h —>■ s'-ff decay via a new class of 
Feynman diagrams in which the W boson is replaced by the charged scalar. However, the 
only important contributions are those with a top quark in the penguin loop, because of the 
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proportionality of the scalar-fermion coupling to the fermion massQ. Therefore, the amplitude 
for b — > S77 in Model I and Model II will be given by Eq. (^) where we only modify T^^ip^i to 



include the new form factor |]18 



36(1 - 

6 - yt) [3 - 8yt + byl + 2(2 - ?>yt) \ogyt] ] , 



i^l™(2/t) = ^4 [7 - I2y, - 3y^ + 8yf + 6y,{2 - 3y,) logy,] + (15) 

36(1 - ytP 



where yt = {mt/MhY and denotes the mass of the charged scalar. Moreover we have used 
the notation^: ^ = fi/f2 = 1/ tan/3 and respectively 



i' = i in Model I (16) 
i' = ~\ in Model II . 

1.2 
1.1 

1 
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Mh (GeV) 

Figure 5: The ratio -Rafb ^ function of the charged scalar mass Mh, for different values of 
tan/3 in Model I (solid: tan/3 = 0.5; long-dashed: tan/3 = 1; short-dashed: tan/3 = 10) and 
Model II (dotted: tan/3 = 0.5; long-dash-dot: tan/3 = 1; short-dash-dot: tan/3 = 10). The 
curves correspond to the cut parameter c = 0.01. 

Summary of comparisons is given in Table ^ and in Figs. ^-|^. In these two figures the differences 

between the 2HDM's and the SM are parametrized in terms of the two ratios 

*This means that the new contribution to the rric and 77ii,-dependent terms in the amphtude are much smaUer 
than the corresponding SM contributions, both IPR and IPI. 

^Note that the notation of Ref. for Model I and Model II is the reverse of ours. 
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Br{B X,77)SM 



and -Ra] 



-PB 




2HDM 
FB 



(17) 



both for Model I and Model II. For illustrative purposes we consider three values of tan/3, 
namely tan (3 = 0.5, 1, 10 and we allow to range between 100 GeV and 1 TeV. In this regard 
we should recall that is also somewhat constrained by the experimental measurement of 
Br{B Xs'y) (see Eq. (|l|)). The present situation seems to indicate that Mh > 200 — 300 
GeV. However both experiment and theory still have an appreciable error. Therefore we are 
tentatively considering > 100 GeV. The upper bound of 1 TeV is, of course, dictated by the 
requirement of a weak-coupled scalar sector. For Figs. § and ^ we also fixed c = 0.01 and kept 
the other cuts on and on the angles as before. In fact a different value of the cut parameter 
(e.g. c = 0.02) does not change the curves in Figs. ^ and ^ significantly. 

We can see that Model II gives always both a bigger branching ratio and a bigger asymmetry 
than the SM. However while the asymmetry is very close to the SM value, the branching ratio 
can be much larger than the SM one over a wide range of values of the scalar mass, especially 
for light Mh. 

The behavior of Model I is quite different. Both the branching ratioQ and the asymmetry 
may be bigger or smaller than the SM, over different ranges of and for different values of 
tan/3. Moreover the relative importance of IPR and IPI contributions may be very different 
than in the SM, due to a cancellation between F2^{xt) and F2^^^{yt) for some particular values 
of Mh and tan (3. As we can see from Table |l], there are cases in which the IPI contribution 
(i.e. mainly the light quark one) is dominant. It is also interesting to note how different values 
of tan (3 imply a totally different behavior both of Rbt and of -Rafb and they complement each 
other by giving substantial deviations from the SM in different regions of the Mh spectrum. 
For example, in the case of tan/3 = 0.5, for light Mh the branching ratio gives a much larger 
deviation from the SM than the asymmetry, while for heavier Mh, even up to Mh — 1 TeV, 
the asymmetry can be more interesting, especially because the branching ratio is smaller than 
the SM one in this region. For tan/? = 1, on the other hand, both the branching ratio and the 
asymmetry are quite different from the SM for Mh < 500 GeV. 

This description has one obvious shortcoming, namely that QCD corrections are not in- 
cluded. However we do not think that the inclusion of QCD corrections in either the SM or 
in the 2HDM calculation should greatly change the important features that we have decribed 
above. 

similar behavior was found in Ref. [p9l for the 6 ^ 57 decay. 
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4 Experimental Considerations 



Our primary focus of course has been on the inclusive b 577 process. As always it is 
very difficult to make reliable predictions about exclusive channels. Nevertheless, some general 
remarks can be made. Unlike 6 — > 57, for b —>■ 577 the pseudoscalar {K for B, rj and rj' for 
Bg) as well as the vector {K* for B and for Bg) final states are both allowed. These single 
meson states are likely to be a large fraction, perhaps several tens of percents of the total 
B Xs77 sample. Since the inclusive branching fraction of charmless B events is roughly 
estimated to be about 1%, they are likely to provide the most important background. However 
it is useful to note that a very important characteristic of 6 — S77 is that the photons carry 
off an appreciable fraction of the total energy, leaving the mean energy of the s quark to 
be around 1.7 GeV. Thus we expect these 27 events to be relatively clean with an average 
multiplicity substantially less than in charmless B events. This property should come in handy 
for separating the 27 background coming from the decay of n^, rj and rj' from amongst the 
multibody charmless B events. Another remarkable feature of the prompt b — > 577 signal is 
that the photons have a large opening angle ((cos^^^) ~ —0.7), as is seen in Fig. ^. This is in 
sharp contrast to the 27's from decays of the 7r°, rj or rj', which, as a rule, should have a much 
smaller opening angle. 
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